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Abstract 

Asymptotic Pade-approximant methods are utilized to estimate the 0{al) contribution to the H gg rate 
and the 0{a1) contribution to the H ^ bb rate. The former process is of particular interest because of the slow 
convergence evident from the three known terms of its QCD series, which begins with an 0(ag) leading-order 
term. The 0{a^) contribution to the H gg rate is expressed as a degree-3 polynomial in L = ln{fM^ / rrii (fi)) . 
We find that asymptotic Pade-approximant predictions for the coefficients of L, , and are respectively 
within 1%, 2%, and 7% of true values extracted via renormalization-group methods. Upon including the full set 
of next-order coefficients, the H —* gg rate is found to be virtually scale-independent over the 0.3 Mh ~ ~ Mt 
range of the renormalization scale-parameter /i. We conclude by discussing the small C'(as) contribution to the 
H ^ bb rate, which is obtained from a prior asymptotic Pade-approximant estimate of the 0{ai) contribution 
to the quark-antiquark scalar-current correlation function. 

1. Introduction 

The Higgs boson characterizes the electroweak symmetry breaking underlying the Standard Model. Its discovery 
and phenomenology will be of immense importance in clarifying our understanding of this symmetry breaking, as 
well as in providing vital information as to the nature of beyond-the-Standard-Model physics. The two leading 
hadronic decay modes of a Weinberg-Salam Higgs boson (iJ) with mass between 100 GeV and 160 GeV are 
the QCD processes H ^ bb and H two gluons (gg). Although the rate for this latter process is known 
to 0{af), such precision incorporates only two non-leading orders of a slowly converging series in the strong 
coupling. If Mh = 100 GeV, for example, the known order-by-order QCD corrections have been calculated ^ to 
be 1-h 0.66 -I- 0.21. 

The work presented here is primarily directed toward obtaining a more precise estimate of the H — > gg decay 
rate. We utilize renormalization-group (RG) and asymptotic Pade-approximant methods to estimate the full next- 
order contribution to the underlying correlation function for this process. Such an approach has already been 
applied to the quark-antiquark scalar-current correlation function underlying the H ^ bb rate , a calculation we 
review in Section 5 of the present paper. 

As in our prior analysis of the two-gluon decay amplitude of a non-Standard-Model CP-odd Higgs field |^ , the 
approach we take here is to test asymptotic Pade-approximant estimates against RG-accessible coefficients within 
the next-order of perturbation theory. We show such estimates to be accurate up to relative errors of only a few 
percent, supporting the credibility of the same approach in estimating the HG-inaccessible coefficient needed to 
determine the full C(q!^) contribution to the H ^ gg rate. 



We operate within the context of MS expressions for the H ^ gg and H ^ bb rates that expHcitly depend 
on an arbitrary renormahzation scale /x. It has been argued elsewhere |^, ^ that asymptotic Pade-approximant 
methods reduce the explicit scale-dependence of perturbative quantities which must ultimately be scale invariant. 
We find this to be the case for the H ^ gg rate as well, despite residual scale-sensitivity anticipated from the 
estimated RG-inaccessible coefficient of the nonlogarithmic 0{a^g) term within the next-order correlator. 

In Section 2, we demonstrate the explicit RG-invariance of the H ^ gg rate, as calculated in in the 
rrib 0, Af|f << 4:Mi limit. This RG invariance enables one to calculate all the next-order coefficients Ck of 
ln'''[/i^/mf (/i^)] (as(/x)/7r)^ within the calculated rate; only the k = constant term is RG-inaccessible. 

In Section 3, we review how asymptotic Pade-approximant methods may be utilized to estimate this set of 
next-order coeflicients, and demonstrate close agreement with the RG-determinations of Section 2 over the range 
100 < Mh < 175 GeV. The RG-inaccessible coeflticient cq is also estimated over this range of Higgs masses and is 
fitted to its anticipated behaviour as a degree-3 polynomial in ln[{Mff/Mi)poie]- 

In Section 4 we examine the residual scale dependence of the Pade-improved H —>■ gg decay rate, as obtained 
in Section 3. Although we anticipate a relative scale dependence comparable to co(as(A*)/7r)'^(~ 3 — 10%), since 
Co cannot be extracted perturbatively from the RG equation (to the order we consider), we find that the residual 
scale dependence of the decay rate is a full order of magnitude smaller than this estimate for /i between 0.3 Mh 
and Mt (we assume 100 GeV < Mh < 175 GeV). Consequently, the scale-dependence of the Pade-derived term 
co(as(/i)/7r)^ very nearly cancels the scale-dependence anticipated from truncation of the perturbative series for 
the H ^ gg rate, thereby facilitating the (nearly) scale-invariant rate predictions of Table 6. Predictions of the 
H gg decay rate are tabulated both for the rrib 0, << AM^ limit, and for the leading-order departures 
from this limit for appropriate constant and running values for mi, and Mt. 

In the concluding section, we review how prior asymptotic Pade-approximant estimates of the 0{a^) contribution 
to the scalar-current correlation function ||^ can be incorporated into the H bb decay rate. We demonstrate 
that the resulting next-order contribution, though only 0.01% of the leading term, is somewhat larger than known 
0{al) power-suppressed contributions |Q. Consequently, the estimated ©(a^) term enables the H ^ bb rate to be 
estimated to four significant figure accuracy. 

2. RG-Invariance of the H ^ gg Rate 

The Higgs gg decay rate is given explicitly to 3-loop order in |Q] by the following expressions: 

Tn^gg = y^R{a,..q' - M^,A^^M2), (2.1a) 

IUh 

R^Cflm<[0[]^ >, (2.1b) 
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In the above expression, < [O'^J^ > is the vacuum polarization of the Higgs field induced via the gluon operator 
§, and = /tt, where a^J^^\ii) is the running strong coupling with n/ active flavours. 
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Five flavours are assumed to be light in both (2.1c) and (2. Id). The i-quark mass Mt appearing in (2.1c) is an 
RG-invariant pole mass, and Mfj is assumed to be small compared to ^M^ [|l|. 

Our goal here is to use asymptotic Pade-approximant methods in conjunction with the RG-invariance of (2.1a) 
in order to predict the next order contribution to T{H — s- gg). In a previous application |^ of such methods to the 
inclusive semileptonic h ^ u rate, it was found that success in predicting RG-accessible next-order coefficients was 
greatly enhanced by recasting the entire expression in terms of the running fermion mass. Indeed, such replacement 
of the 6-quark pole mass with its scale-dependent MS mass had already been employed by van Ritbergen |^ to 
avoid a renormalon pole. 

Consequently, we first re-express the H 2g rate in terms of the running t-quark mass mt{fj,), which evolves 
via a six-active- flavour 7r„-function, and the corresponding six-active-flavour running coupling x^^\ii). This trans- 
formation is facilitated by the following relationships M : 



(m) 



6 



Ir 



O 



(2.2) 



mt{fi)/Mt = 



r.(6)/ 



(^) ( O + ^" ( 



O 



(a;(6))2 



(2.3) 



Equation (2.3) implies that the logarithm in (2.1c) may be re-expressed in terms of L = Zn(/i^/mj (/i)): 



( T72 ) = ^ - ^^^^Ha^) [2i + 8/3] + O (a:(«0 



(2.4) 



With = Mfj, the logarithm in (2. Id) can also be expressed entirely in terms of the running i-quark mass and a 
logarithm T = In {^M^/Mf^ of the ratio of iJG-invariant pole masses: 



Ml 

Substitution of (2.2) and (2.5) into (2. Id) with 
expression for the H gg decay rate: 
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Mfj, and substitution of (2.4) into (2.1c) leads to the following 
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In (2.7), we list the unknown coefficients co, ci, C2 and C3 of the 4-loop contribution to the rate. Three of these 
may be extracted by the scale-[RG-] invariance of the physical decay rate: dT /dji = 0. This invariance implies that 



dS 

= [l-27™(.%+/3(x)^. 
Both the (3 and 7,„ functions in (2.8) are referenced to six active fiavours: 



(2.8) 
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ji^){x) = -x--x'... . (2.10) 

One can easily verify from the known terms listed in (2.7) that (2.8) is perturbatively valid to orders , (including 
explicit cancellation of terms involving T), and x'^L. The continued perturbative validity of (2.8) to orders x^L^, 
x^L, and x^ is sufficient to determine the four-loop coefficients ci,C2, and C3: 

910.3167 -i|15T + 5I25^^ (2.11, 

1225 1127 

C3^f. (2.13) 

The coefficient cq is not RG-accessible to these orders. In the next section, we will utilize asymptotic Pade 
approximant methods to estimate the four-loop coefficients {cq, ci, C2, C3}. As in prior work ^, the accuracy 
of these predictions in reproducing (2.11 - 2.13) will serve as an indication of the accuracy of our estimate for cq. 

3. Pade-Predictions for the Four-Loop 
Coefficients 

The asymptotic Pade-approximant procedure for estimation of the four loop coefRcients cq, 01,02,03 in (2.7) has 
been delineated in previous work ||, |^. The series (2.7) may be expressed in the form 

S[x,L,T]^x^ [1 + Ri[L,T]x + R2[L,T]x^ + R3[L,T]x^ + ...] (3.1) 

where i?i and R2 are known, and R3 is to be determined: 

^ ' ' 48 48 V 16 8 y 16 ^ ' 

RsiL, T] = co(T) + oi{T)L + 02^1^ + cgL^. (3.4) 

Initially, we shall eliminate T as a variable by assuming that the Higgs pole mass is 100 GeV, in which case 
T = InlMjj/Mf] = 2 /ri(100/175.6) = -1.126. As described in §] and the [0|1] Pade-approximant prediction 
for i?2 is 

4""' = R\ (3.5) 
and the [1|1] Pade-approximant prediction for i?3 is 

4^l'l = Rl/R^. (3.6) 

If the error of [N\l] approximants in predicting RN+2^ the N + 2 term in the perturbative series, is inversely 
proportional to -I- 1 liol — i.e., if 
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where ^ is a constant — one may then utilize (3.5) and the exact value for i?2 within (3.7) to obtain A — {R2 — 
Ri)/R2- Substituting (3.6) and this estimate for A into (3.7), one obtains an error-improved estimate for the 
unknown coefhcient R3 

R3[L] - 2Rl[L]/{Ri[L]R2[L] + RI[L]) (3.8) 

To obtain estimates of the coefficients co,ci,C2,C3 within (3.4), we match the scale dependence of (3.4) to that of 
(3.8) over the purely perturbative L > region [corresponding to the ultraviolet scales /x > TOt(/x)] through use of 
the moment integrals 

Nk = {k + 2) f dww''+^R3{w), (3.9) 
Jo 

where w = mf{fi)/fi'^ [L — ~ln{w)]. Substitution of (3.4) into the integrand of (3.9) yields the following expressions 
for the first four moments 

N_i ^ Co + ci + 2c2 + 6C3, (3.10a) 
113 

No = co + -ci + -C2 + -C3, (3.10b) 
12 2 

Ni^co + -ci + -C2 + -C3, (3.10c) 

^2 = Co + + + ^C3. (3.10d) 

However, explicit numerical estimates of these four moments may be obtained via substitution of (3.2) and (3.3) 
for the respective factors of and R2[L] appearing in (3.8), and by subsequent substitution of this estimate 

for R^lL] into the integrand of (3.9) [with L = —ln{w)]. For Mjj = 100 GeV the resulting estimates are 

A^_i = 5102.9, iVo = 3542.9, Ni ^ 3131.8, N2 = 2944.7. (3.11) 

Substitution of these values into (3.10) yields the following predicted values {cf'^'^^) for the four- loop terms in the 
H ^ gg decay rate (2.7): 

^Pade ^ 2453, c^'"^^ = 1772, c^"'^^ = 378.4, c^"'^^ = 20.27. (3.12) 



Cn 



The exact values of ci,C2, and C3 were determined via RG methods in the previous section. For Mh =^ 100 GeV 
[i.e. for T = —1.126], these values are found from (2.11-13) to be 

ci = 1789.02, C2 = 385.568, C3 = 21.4375. (3.13) 

Comparing (3.12) and (3.13), one finds the relative error [Sci = {cf'^'^'^ — Ci)/ci\ of the Fade estimates for ci, C2 and 
C3 to be -0.95%, -1.9%, and -5.5%, respectively. 

Such accuracy in predicting the known four-loop terms in the H gg rate suggests that the prediction for cq in 
(3.12) is a credible one. One way to test the stability of this prediction is to utilize the true values of ci, C2, C3 within 
the moment expressions (3.10) to obtain [for the numerical values (3.11) already obtained for these moments] four 
independent determinations of cq. We then find that 
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results all within 2% of that in (3.12). 

An alternative approach to matching the four-loop coefficients within (3.4) to the error improved Pade 
estimate (3.8) is optimize the least-squares function ||^ 

X'^[co, ci, C2, C3] = / [R3 — {cq — cilnw + C2ln^w — c^ln'^w)]'^ dw, (3.15) 



with i?3 in the integral given by (3.8). As before, factors of Ri and R2 appearing in (3.8) are given explicitly by 
(3.2) and (3.3) with L = -ln{w). One then finds for Mh = 100 GeV [T = -1.126] that 

X^(co, ci, C2, C3) = 4.064168878 • 10^ + 720c^ + 4coC2 + + 2coCi + 24c^ + 12ciC2 + 12coC3 + 2c? 

-I- 240C2C3 + 48ciC3 - 10205. 87co - 17507.05ci - 54106. 23c2 - 234536.67c3 (3.16) 

The optimization requirement 

(3.17) 

yields predictions remarkably close to those of (3.12), 

cf = 2452, cf = 1774, ^ 377.2, cf = 20.45, (3.18) 

further confirming the stability of the estimation procedure. 

In Table 1, we have tabulated a set of predictions of the four-loop term cq for values of the Higgs mass between 
100 and 175 GeV. Also tabulated are the errors in the predicted values for ci, 02,03, relative to the true values 
for these coefficients, as given in (2.11-13). Estimated values of ci and C2 remain within 2% of their true values 
(2.11-13) over the range of Higgs masses given; the relative error of C3 estimates remains below 7% for the same 
range. This consistency provides further support for the cq estimates presented in the final column. 

These cq estimates may be utilized to ascertain the convergence of the series (2.7). If we choose ^ — mt{fj.), 
the logarithmic factors within (2.7) all vanish. To evaluate the series S/x'^ in (2.7), we assume that mt{mt) = 
Mt = 175.6 GeV and that as(175.6) = 0.10915 = 7ra;(175.6), as evolved from a^(M^) = 0.119 |l|, in which case 

S{fi = 175.6 GeV)/a;2(175.6) = 1 + aox(175.6) + box^{175.6) + 00X^(175.6). (3.19) 

For Higgs masses between 100-175 GeV, we tabulate in Table 2 the magnitudes of successive orders in (3.19). The 
final 4*'* order term is obtained from the appropriate Pade estimate for cq listed in Table 1. Table 2 shows that 
the four-loop term decreases from 10% to 3% of the leading contribution as Mh increases from 100 to 175 GeV. 
Moreover, the convergence of the series is problematical in the absence of the estimate for the four loop term. Over 
the range of Higgs masses considered, the three-loop contribution is between 18% and 33% of the leading (one-loop) 
contribution. 

The procedures delineated above can also be utilized to predict cq's explicit polynomial dependence of T, 

co{T) = ao + aiT + a2T'^ + aaT^ (3.20) 

analogous to the expressions (2.11-13) obtained from _RG-invariance for ci{T), C2{T), and C3(T). To extract the T 
dependence of cq, we first incorporate all known T-dependence into the least- squares function (3.15): 



x'[co(T)]=/ [R3[L,T]-{co{T) + ci{T)L + C2{T)L^+C3{T)L''}Ydw. (3.21) 
Jo 

In (3.21), L = —Inw. The quantities ci(T), C2(T) and C3(T) are no longer optimizable variables as in (3.15), but 
are now the explicit polynomials (2.11-13) obtained in the previous section via i?G'-methods. The factor R^lL^T] 
in the integrand of (3.21) is just (3.8) generalized to include the explicit T-dependence of Ri and R2 [eqs. (3.2) 
and (3.3)]: 

^1 = R, [L, TmtrH Rl [L,Ty ^^-^^^ 
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From (3.21), the requirement dx^ /dco — generates cq as a function of T. Since T = ln{M'^ / M^) , we restrict our 
attention to the region —1 < T < (i.e. to values for Mh between Mt and Alte^^^^ = 107 GeV). A set of values 
for co{T) can be obtained via optimization of (3.21) over values of T in this region: 

co(0) = 735.7, co(-O.l) = 841.7, co(-0.2) = 955.6, co(-0.3) = 1078, 
co(-0.4) = 1208, co(-0.5) = 1346, co(-0.6) = 1493, co(-0.7) = 1649, 

co(~0.8) = 1814,; co(-0.9) = 1988, co(-l.O) = 2170. (3.23) 
We obtain a least-squares fit of these results to the form (3.20) by optimizing 

10 

X^[ao,ai,a2,a3] = ^ 

1=0 

with respect to {oq, ai, 02, 03}, and find that 

co(T) = 735.7 - 1020T + 388.8T^ - 25.4ir^ (3.25) 

Each coefficient listed above is within 5% relative error of the corresponding coefRcient obtained via a least squares 
fit of the Co values displayed in Table 1 to the polynomial form (3.20): 

co(T) = 755.9 - 1029r + 3M.3T^ - 26.7AT^. (3.26) 

Note that (2.11), (2.12), (2.13) and the Pade-estimates (3.25) or (3.26) specify all the logarithmic coefficients 
within the full four loop contribution (3.4) to the H ^ gg decay rate (2.6,7), and comparison of these results 
with future perturbative calculations should yield information which can be employed to further improve Pade 
estimation procedures. The exact H ^ gg rate (2.6) is necessarily a scale invariant physical quantity. The factor 
S[x{^), L{p),T] within (2.6) will exhibit residual scale-dependence [i.e. /^-dependence] only as a consequence of 
truncation of the series S* to a given order of perturbation theory. In the section that follows, the incorporation of 
the four-loop coefficients co_3 of S[x{p)^ L{ii), T] is seen to eliminate virtually all of this residual scale dependence. 
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4. The H ^ gg Rate 

In Fig. 1, the /i-dependence of three- and four-loop expressions for S[x''^\^), i(/i),T] is plotted for the case of 
a 140 GeV Higgs mass [T = 2 Zn(140/175.6)]. The four-loop term within S is evaluated through use of the Pade 
estimate cq = 1306 (Table 1) in conjunction with eqs. (2.11-2.13) for the iJG-accessible coefficients ci,C2, and 
C3. The running coupling as{^) and running mass mt{ij) occurring within (2.7) are evolved via four-loop (3 and 
7-functions from physical reference values as(Mz) = 0.119 and mt{mt) — 175.6 GeV jl^ . 

Figure 1 shows that the four- loop expression eliminates virtually all of the residual scale dependence still evident 
in the three-loop rate in the region > /i > 30 GeV. The four-loop rate is observed to have a local minimum at 
fj, — 43.5 GeV and a weak local maximum at = 89.5 GeV. The values of S[x(p), L{p), T] at both of these points 
of minimal-sensitivity differ by only 0.2%, indicative of the flatness of S between these two points. In Table 3, 
values for S[x{ii), L{fi), T] as well as the term-by-term series = x^[l + Rix + R2X^ + Rsx^] within S are displayed 
for a variety of /(-values of interest. The table displays the relative size of successive terms at the minimal- 

sensitivity points /i — 43.5 and 89.5 GeV in addition to the points /i — 140 GeV(= Mh) and n = 175 GeV(= Mt). 
The relative magnitude of the four- loop term Rsx^ is seen to be less than 1% of the leading term in the series 
(unity) over the entire region between the minimal-sensitivity values of fi. 

Of particular interest in this range are those values of fi [47.0 GeV and 73.5 GeV] at which the four-loop term 
R^x^ effectively vanishes. These values correspond to the two points in Figure 1 at which the 3-loop and 4-loop 
curves cross. The usual approach towards extracting information from an asymptotic series X]n=o ^nx"' {Rq — 1) 
is to sum only that series' decreasing terms, i.e. to evaluate X]n=o choosing n' such that \Rn'X^ \ is 

a minimum. By choosing ^ so as to have R^^fJ-) vanish, one can then argue (for this choice of /() that n' — 3, 
suggesting that such a value for ^ is optimal for estimating the series from its known terms. Of course, such an 
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interpretation rests on the assumed increase of terms \Rnx"-\ subsequent to n' = 3; all we can really be certain of 
is that |i?4a:''| > |i?3a;^| (= 0) at such a value of fi. Nevertheless, Table 3 shows that both values of /i for which 
R3 —^ lie between the two minimal-sensitivity points [/i = 43.5 and 89.5 GeV], and that the /i-sensitive factor 
S[x{fi), L{ii), T] within the rate (2.6) varies by less than 0.2% over this entire region. Over the full range of ^ values 
displayed in Table 3 (43.5 GeV < fi < 175 GeV), the four-loop series term (R^x^) varies between and 5% of the 
leading one-loop order term (unity). Surprisingly, however, the rate S[x{^), L{fj,),T] displays a relative spread of 
values AS/ S < 0.4% over the same range of fj,, indicative of a substantial reduction in residual scale dependence. 

In Figures 2 and 3 we exhibit the residual scale dependence of S[x{fj,), L{^),T] to three and four- loop order for 
Higgs masses of 100 GeV and 175 GeV, respectively. These figures show the same reduction in scale dependence 
evident in Figure 1 when Mh = 140 GeV. For the case of Mh ~ 100 GeV, for example, Table 4 shows that the 
4-loop term R^x^ varies between zero and 10% of the leading-order 1-loop term between the local minimum at 
/Lt = 29 GeV and /i = Mt. Over this same range of /i, the four-loop values of S[x{^), L(/x), T] remain within 1% of 
each other. For = 175 GeV, the four-loop term is between zero and 3% of the leading-order term between the 
local minimum at /i = 48.5 GeV and fj, = Mt, whereas the full four- loop expression for S[x{fi), L{p), T] exhibits a 
relative spread of only 0.1% [Table 5]. 

Such scale independence is, of course, a reflection of the i?G'-invariance (2.8) of S[x{^), L{^),T], which has 
been utilized explicitly to obtain the coefficients ci_3. Nevertheless, the coefficient cq, which is not perturbatively 
accessible via (2.8) but is obtainable (at present) only by Fade approximant methods, appears to be precisely what 
is required to eliminate virtually all residual /i-dependence in the rate arising from truncation. For example, if 
Mh = 100 GeV, the factor cqx^{^) in isolation contributes more than 10% of the leading order contribution (unity) 
for J; Mt- Moreover, this contribution increases as /x decreases. Nevertheless, the particular choice cq = 2453 
appears to ensure that the overall spread in 5* remains within 1%, despite the potentially large contribution to this 
spread from cox'^(/i). Evidently the /x dependence of this term serves to cancel the residual /i-dependence of the 
remaining terms in the series. |^ 

This near cancellation of residual scale-dependence makes possible a set of credible predictions for the rate 
(2.6). In the fourth column of Table 6 we have tabulated the H ^ gg decay rates for Mh = 100, 125, 140, 150, and 
175 GeV in the mb 0, Mh << 4:Mf limit. The largest source of uncertainty for these predictions is in the value 
for as{Mz) — (0.119 ± 0.002 ||l2|), which should lead to 4% uncertainty in the rates presented in Table 6. 

Fermion mass effects {i.e. the departures from the mi, — )■ 0, Al'^ << AM^ assumptions implicit in the derivation 
0] of the H gg rate) can be accommodated in leading order by replacing the factor of unity in (3.1) with the 
following m^- and Afj-sensitive terms |l^ : 
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(4.4) 



The right hand side of (4.1) is easily seen to approach unity when mb 0, Mfj/AM^ — > 0: 



lim At — —] lim Ab = 0. 

Tt^O 3 Tb^OC 



(4.5) 



The leading mass correction (5,„, as defined in (4.1), is tabulated for various Higgs boson masses (with PDG [|T2j 
fermion-mass values mb = 4.2 GeV, Mt = 175.6 GeV) in the fifth column of Table 6. This lowest-order fermion 
mass correction can be incorporated into the rate (2.6) by its inclusion into the series (3.1) 



5 = a;2[H- Rix + R2X^ + R^x^ + ^™]. 

^ The reduction of scale-dependence via Pade approximant methods is discussed in detail by Gardi 



(4.6) 
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The sixth column of Table 6 tabulates the H ^ gg rate with this correction included. 

Because 5m is a leading order correction, there is genuine ambiguity as to whether "physical" or running masses 
should be incorporated into this correction. The former choice yields a manifestly scale-dependent contribution 
AS* = SP^x'^{^) to the rate S. One could argue for the incorporation of running masses TOf,(/^), mt(/i) in (4.2) and 
(4.3); i.e., AS = (5m(/x)a;^(/i), with (5m(/x) calculated via (4.1-4) with Tf = Mfj/im^di) and Tf, = Mfj / imKii) , and 
with /i identified consistently with the minimal-sensitivity scale tabulated in Column 2 of Table 6. In Column 7 
of Table 6 we have tabulated SmifJ-) utilizing 6-active-flavour running masses referenced to mt{mt) — 175.6 CeV, 
mb{mb) = 4.2 CeV, and in Column 8 we have listed the corresponding H — > gg decay rates. Theoretical uncertainty 
associated with the masses utilized in the leading correction-factor 5m is reflected in the differing rates of Column 
6 and Column 8. This discrepancy is seen to be at most a 2% effect. 

The leading order mass-correction factor 5m in (4.6) is seen from a comparison of Table 6 to Tables 3-5 to 
be generally smaller than the three-loop term R2X^ in (4.6), but somewhat larger than the four-loop term R^x^ 
also appearing in (4.6) [5m is smaller than R-iX^ when /i = Mt, as in Table 2]. Since the next order of fermion 
mass corrections is suppressed by an additional power of as , we anticipate such next-to-leading-order fermion mass 
corrections to be well within 1% of the total rate. In any case, both leading-order mass corrections 5m and the 
four-loop contribution R^x^ we obtain are sufficiently small for credible estimates of the H ^ gg rate on the 
basis of the terms in (4.6). It is also worth noting from Table 6 that 5f^ itself is only 1% of the leading-order one 
loop term [normalized to unity in (4.6)] when Mh = 140 — 150 GeV, and that 5m{lJ-) is comparably small when 
Mh = 150 CeV. Hence, the rates tabulated in Table 6 are optimally valid for Higgs masses in the 140-150 GeV 
range. 



5. The H ^bb rate 



The leading QCD scalar-current correlation function contributions to the Higgs — > bb decay rate are known from 
explicit calculation to ©(a^), with secondary 0{ni\/ Mfj) power corrections known to 0{a1), as given in 



r(iJ ^ bb) 



4\/27r 



MnmliMH) 



^scalar {fJ-^ = S = Alfj) 

mliMn) 



+- 



Ml 



_6- 40^^^^-87.72 



as{MH) 



(5.1) 



The factor liscaiar in (5.1) is the imaginary part of the QCD correlation function for the quark-antiquark scalar 
current normalized to unity in leading order Q: 



nsca/ar[M- S,x{^i)] = 1 + 

+ 
+ 



17 

y + 2 In^i^/s ] x{fj,) 



263, 



47,^2 ^A^' 



29.1467+ — + — ^ 
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41.7576+ 238.381 In ( ^ 



94.6759 In^ ( — 



7.61574 Zn^ [ ^ 



x^p) 



diln ( — 



d^ln^l^] ) x\n) 



(5.2) 



where a:(/z) = [as{n)]nj=5/''T, corresponding to five active flavours. 

In ref.[2], a detailed asymptotic Pade-approximant procedure is presented for estimating the coefhcients do-i 
in (5.2). The methodology is virtually identical to that of Section 3 above, except that five moments are now 
calculated for R4 [w] , the coefficient of x'^ (fj,) , based upon the asymptotic error formula prediction ||Tl[| 

+ R1R2R3 — 2i?i^i?3] 



i?4 — 



(5.3) 
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The results, as tabulated in Table 3 of ref. [2], are 

do 64.2, di = 745, d2 = 1180, = 253, di = 15.4. (5.4) 
The factors c?i_4 may be obtained directly from RG-invariance of the physical rate 

[mli^ijUscaiar [^i, s, a;(/i)]] = 0, (5.5) 

leading to the following values (also tabulated) in |^): 

di = 791.52, d2 = 1114.7, d^ = 260.06, d^ = 14.755. (5.6) 

The strong agreement between (5.6) and (5.4) suggests that the estimate for do in (5.4) is a credible one. It is 
interesting to note that a much more naive Pade approach for estimating d^, in which (5.3) was applied directly to 
^scalar[^J'^ = s = Mfj), yielded a value for do {—67.25 [|ll]) surprisingly consistent with the estimate in (5.4). Only 
the latter estimate exhibits sensitivity to the logarithmic terms in Tlgcaian which all vanish when fi^ — s. 

If we utilize the estimate for dg quoted in (5.4) within the H ^ bb rate, we find that the correlation-function 
factor within (5.1) is given to ©(af ) by 

n,ca/ar(Ai' = s = Aljj) = 1 + ya;(AfH) + 29. 1467 (Mjf ) + 41 .7576 a;3(Afff ) + 64 a;^Afo). (5.7) 

The underlined coefficient is, of course, do as estimated in ||^ via asymptotic Pade-approximant methods. 

There is genuine value in having an estimate of this term, as it is generally larger than the 0{a^ml/M'^) final 
term in (5.1). For example, suppose that Mh = 130 GeV, mf,(130 GeV) = 2.7 GeV and x(130 GeV) = 0.114/7r 
(these values are consistent with those employed in Table 1 of ref. 0). The relative magnitudes of the correlator 
contributions 95.5) to the H ^ bb rate are seen to be 

Hscaiar (a^^ = s = (130 GcV)^) = 1 + 0.2056 + 0.0384 + 0.0020 + 0.00011 . (5.8) 

All but the underlined Pade-estimate term are tabulated in Table 1 of ref. . Comparison to the power-suppressed 
terms in (5.1), 

^''^"^ (-6 - 40 x{Mh) - 87.72 x^{Mh)) — * - 0.0026- 0.00062- 0.00005, (5.9) 

^ J\/fj = 130GcV 

(also tabulated in 0|) reveals that the 0{a'^) term in (5.8) is double the magnitude of the 0{a^) term in (5.9). 
Hence the Pade-estimated 0{af) term in (5.7) enables one to utilize the full precision available in the known 
power-suppressed contributions to (5.1). However, one cannot anticipate such precision experimentally for many 
years to come. 



6. Summary 

In the preceding sections, we have demonstrated how asymptotic Pade-approximant methods may be utilized to 
estimate the unknown O (a^) coefficients {cq, Ci, C2, C3} within the H ^ 2g decay rate (2.7). Such estimates for 
{ci, C2, C3} are seen (Table 1) to be within a few percent of the true values for these coefficients, which can be 
extracted via renormalization-group methods. This accuracy supports corresponding asymptotic Pade-approximant 
estimates of the renormalization-group- inaccessible coefficients co presented in Table 1. Moreover, the inclusion 
of O (og) terms within (2.7) is seen to virtually eliminate the scale-parameter dependence of the rate over an 
astonishingly large range of the scale-parameter /i, typically O.SMh i fi ^ Mt- Inclusion of estimated O (a^) 
corrections in conjunction with leading-order fermion-mass corrections to the rate (Table 6) are seen to reduce the 
perturbative uncertainty in the H 2g decay rate from 0(20%) to 0(2%). Four-loop corrections to the H ^ bb 
decay mode are also presented, which are seen to reduce the very small perturbative uncertainty of this dominant 
hadronic mode by an additional order of magnitude. 
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Mjj(GeV) 








Co 


100 


-0.010 


-0.019 


-0.054 


2453 


125 


-0.009 


-0.013 


-0.059 


1646 


140 


-0.009 


-0.009 


-0.062 


1306 


150 


-0.009 


-0.007 


-0.068 


1120 


175 


-0.012 


-0.0004 


-0.070 


763 



Table 1: The final column tabulates asymptotic Pade-approximant predictions for the four loop nonlogarithmic 
coefficient cq within (2.7), as obtained via moments of cq — cilnw + C2ln'^w — c^ln'^u!, for Higgs masses between 
100 GeV and Mf. Relative errors between predictions and true values of ci,C2 and C3, as determined from RG 
invariance, are tabulated in the second, third and fourth column: 6ci = ( (J"''=*'=*«='^ _ cf'^ j /cf^. 



Mh(GcV) 


S[x{Mt),L{Mc).T]/x'^ = l + aox + b^x"^ + cqx^ 


100 


1 + 0.7728 + 0.3333 + 0.1014 


125 


1 + 0.7133 + 0.2675 + 0.0679 


140 


1 + 0.6831 + 0.2361 + 0.0538 


150 


1 + 0.6647 + 0.2177 + 0.0461 


175 


1 + 0.6237 + 0.1783 + 0.0312 



Table 2: The perturbative series S/x"^ within the H ^ 2g decay rate (2.7), evaluated at = mt{ij) = 175.6 GeV. 
a;(175.6GeV) = 0.10915/7r, consistent with x{Mz) = as{Mz)/'K = 0.119/7r. The final underlined term of the series 
is obtained from the cq estimates in Table 1. 



/x(GeV) 


a,,(/i) 


mt{p){GcV) 


1 + Rix + R2x'^ + Rsx^ 


S[x{fi),L{fi),T] 


43.5 


0.1325 


198.2 


1 -f- 0.3811 - 0.0395 -1- 0.0023 


0.002390 


47.0 


0.1309 


196.8 


1 -1- 0.4014 - 0.0251 - 0.0001 


0.002390 


73.5 


0.1226 


188.8 


1 + 0.5094 + 0.0614 + 0.0000 


0.002393 


89.5 


0.1193 


185.6 


1 -t- 0.5525 -1- 0.1007 + 0.0068 


0.002394 


140 


0.1124 


178.8 


1 -1- 0.6419 + 0.1907 + 0.0344 


0.002390 


175 


0.1092 


175.6 


1 -1- 0.6825 -1- 0.2354 + 0.0534 


0.002384 



Table 3: The four-loop expression for the scale-dependent factor S[x{iJ,), L{n),T] for Mh = 140 GeV within the 
Higgs gluon-gluon decay rate (2.6) for five different values of fi. Weak extrema within the rate (Fig. 1) occur 
at = 43.5 and 89.5 GeV. Also displayed are /x = Mh, ji Mt, and two values of /x (47.0 GeV and 73.5 GeV) 
at which the four-loop contribution to S is almost zero. The values of S over this entire range of n are seen to be 
remarkably static, as discussed in the text. 



Ai(GeV) 




mtifi) (GeV) 


1 + Rix + R2x'^ + R-iX^ 


S[x{p),L{p),T] 


29.0 


0.1413 


206.5 


1 + 0.3819 - 0.0546 + 0.0033 


0.002691 


31.0 


0.1397 


205.1 


1 -1- 0.4007 - 0.0410 -1- 0.0002 


0.002691 


54.0 


0.1282 


194.2 


1 + 0.5418 + 0.0775 + 0.0001 


0.002698 


65.5 


0.1247 


190.8 


1 + 0.5853 + 0.1198 + 0.0089 


0.002699 


100 


0.1175 


183.8 


1 + 0.6722 -1- 0.2125 -|- 0.0402 


0.002693 


175 


0.1092 


175.6 


1 + 0.7722 -1- 0.3326 -|- 0.1010 


0.002667 



Table 4: The scale-dependence of the four-loop expression for S'[a;(/i), L(/x), T] is displayed, as in Table 3, for the 
case Mh = 100 GeV. Weak extrema in Figure 2 occur at /x = 29.0 and 65.5 GeV. The four-loop contribution is 
seen to nearly vanish at /x = 31.0 and 54.0 GeV. 
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Ai(GeV) 




mtip) (GcV) 


1 + RiX + i?2x2 + RsX^ 


S[x{fi),m,T] 


48.5 


0.1303 


196.2 


1 + 0.3385 - 0.0584 + 0.0081 


0.002216 


61.0 


0.1260 


192.0 


1 + 0.3976 - 0.0190 - 0.0001 


0.002216 


91.0 


0.1190 


185.3 


1 + 0.4912 + 0.0537 + 0.0000 


0.002218 


111 


0.1159 


182.2 


1 + 0.5338 + 0.0910 + 0.0058 


0.002218 


175 


0.1092 


175.6 


1 + 0.6230 + 0.1777 + 0.0309 


0.002215 



Table 5: The scale-dependence of the four-loop expression for S[x{ij), L{ij.),T] is displayed, as in Table 3, for the 
case Mh = 175 GeV. Weak extrema in Figure 3 occur at /x = 48.5 and 111.0 GeV. The four-loop term is seen 
to nearly vanish at /x = 61.0 and 91.0 GeV. The relative spread in the rate S is seen to be only of order 0.1% for 
values of /U between 48.5 and 175 GeV. 



Mh 


M 


S 


^ H 


.=0 




SP 


^ H 


^99 






pi. 

^ H 


-^99 


100 


65.5 


0.00270 


1.97 


10- 


-4 


-0.085 


1.87- 


10- 


-4 


-0.059 


1.90 


10- 


125 


80.5 


0.00249 


3.55 


10- 


-4 


-0.040 


3.47- 


10- 


-4 


-0.015 


3.52 


10- 


140 


89.5 


0.00239 


4.79 


10- 


-4 


-0.012 


4.76 • 


10- 


-4 


-hO.Oll 


4.82 


10- 


150 


95.5 


0.00234 


5.75 


10- 


-4 


-h0.0076 


5.78- 


10- 


-4 


-fO.023 


5.83 


10- 


175 


111 


0.00222 


8.67 


10- 


-4 


-F0.060 


8.99- 


10- 


-4 


-1-0.077 


9.08 


10- 



Table 6: The four-loop order H ^ gg decay rate (2.6) is tabulated for various choices of Mh- The second column 
lists the minimal-sensitivity values of /U (dS/dfi = 0) closest to Mh- The third column is the value of S[x{iJ,), L{^), T] 
evaluated at this choice for /i, as indicated in the previous three tables. The fourth column is the rate in the 
rrib 0, M]j << AMf limit. The fifth column tabulates the magnitudes of the leading fermion mass correction 
to S/x^, as discussed in the text, using "physical" (p) fermion masses Mt = 175.6 GeV, rrib = 4.2 GeV. The sixth 
column tabulates the H gg rate incorporating this mass correction. The seventh column lists the leading fermion 
mass correction to S/x"^ utilizing running fermion mass values mb{iJ,),mt{i^) at the Column 2 values of /x, and the 
final column tabulates the H ^ gg rate incorporating this fermion-mass correction. All masses and decay rates are 
in GeV units. 
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Fig.1 




0.0016 
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0.0012 - 

0001 20 40 60 80 100 120 140 160 

Figure 1: Comparison of the residual renormalization scale- (/i— ) dependence of the three- loop (3L) and four- 
loop (4L) expression for the H ^ gg decay rate of a Higgs boson with mass Mh = 140 GeV. The y-axis numbers 
correspond to the dimensionless scale-dependent factor S [x (/i) , L (/x) , T] within the rate (2.6), and the x-axis is the 
scale fj, in GeV. The logarithmic coefhcicnts of 4L-contributions are extracted via renormalization-group equation 
methods. The 4L non-logarithmic coefficient is obtained via asymptotic Fade approximant methods. The 3L and 
4L curves nearly coincide in the region 44 GeV ^ /i ^ 80 GeV, indicating a very small 4L contribution over this 
range of /x. 
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Fig.2 
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1^ 

Figure 2: A comparison of the /x-dependence of three-loop and four-loop expressions for the H ^ gg rate is plotted, 
as in Figure 1, but now for Mh = 100 GeV. The y-axis numbers correspond to the dimensionless scale-dependent 
factor S [x {ji) ,L{fi), T] within the rate (2.6), and the a;-axis is the scale n in GeV. The curves cross at values of 
fi for which the (estimated) four-loop contribution to the rate is zero. 
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Fig.3 
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Figure 3: A comparison of the /i-dependence of three-loop and four-loop expressions for the H gg rate (as in 
Figs. 1 and 2), for Mh = 175 GeV. The y-axis numbers correspond to the dimensionless scale-dependent factor 
S [x (/Lt) ,L{i^), T] within the rate (2.6), and the x-axis is the scale fx in GeV. 
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